An efficient numerical method for gyrokinetic δ f Vlasov simulations in helical plasmas
نویسندگان
چکیده
Three-dimensional configurations of helical plasmas bring numerical difficulties to gyrokinetic simulations of plasma turbulent transport in toroidal fusion devices. Watanabe et al. have investigated the ion-temperature-gradient (ITG)-driven turbulence in helical plasmas by means of the gyrokinetic Vlasov simulation code GVK [1]. Because of a large number of grid points along a magnetic field line for resolving particle motions in helical ripples, however, the CourantFriedrichs-Lewy (CFL) condition on the parallel advection term as well as the perpendicular magnetic drift term severely restricts time steps in the linear limit. To address this problem, we propose an efficient numerical method for linear analysis of instabilities in helical plasmas by means of gyrokinetic δ f Vlasov simulations, employing semi-Lagrangian and additive semiimplicit Runge-Kutta schemes. The new scheme is free from the CFL restrictions for the linear terms.
منابع مشابه
Particle-in-cell simulation with Vlasov ions and drift kinetic electrons
There are certain limitations in using gyrokinetic ions for simulations of turbulent transport in tokamak plasmas. Applications where Vlasov ions might be more appropriate include the electron temperature gradient driven turbulence, edge turbulence with steep density gradient, and magnetic reconnection in a weak guide field. In such cases the ion gyrokinetic model presently used in simulations ...
متن کاملCollisions in Global Gyrokinetic Simulations of Tokamak Plasmas Using the f Particle-In-Cell Approach: Neoclassical Physics and Turbulent Transport
The present work takes place within the general context of research related to the development of nuclear fusion energy. More specifically, this thesis is mainly a numerical and physical contribution to the understanding of turbulence and associated transport phenomena occuring in tokamak plasmas, the most advanced and promising form of magnetically confined plasmas. The complexity of tokamak p...
متن کاملA Numerical Method for Parallel Particle Motions in Gyrokinetic Vlasov Simulations∗)
A semi-Lagrangian scheme is applied for the first time to computations of charged particle motions along magnetic field lines, to numerically solve the δf gyrokinetic equations in a flux tube geometry. This new solver adopted in the gyrokinetic Vlasov simulations has an advantage over the conventional Eulerian codes in calculating the parallel dynamics, because semi-Lagrangian schemes are free ...
متن کاملGuiding-center simulations on curvilinear meshes
The purpose of this work is to design simulation tools for magnetised plasmas in the ITER project framework. The specific issue we consider is the simulation of turbulent transport in the core of a Tokamak plasma, for which a 5D gyrokinetic model is generally used, where the fast gyromotion of the particles in the strong magnetic field is averaged in order to remove the associated fast time-sca...
متن کاملMagnetic holes in plasmas close to the mirror instability
Non-propagating magnetic hole solutions in anisotropic plasmas near the mirror instability threshold are investigated in numerical simulations of a fluid model that incorporates linear Landau damping and finite Larmor radius corrections calculated in the gyrokinetic approximation. This FLR-Landau fluid model reproduces the subcritical mirror bifurcation recently identified on the Vlasov-Maxwell...
متن کامل